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1 Introduction 

The generalized Funk transform is a integral transform acting from densities on a manifold 
X to functions defined on a family £ of hypersurfaces in X. The dual operator M° is 
again a Funk transform which is defined for densities on the manifold S. We state two- 
side estimates for the operator M and a description of the range of the Funk transform 
operator M and approximation theorem for the kernel of the dual operator. This operator 
is similar to the integral transform related to a double fibration in the sense of Guillemin 
[I] and some results can be extracted from his theory. Other results can be generalized 
for the general double fibration. 



2 Geometry 



Let X and £ be smooth n-manifolds, n > 1 and F be a smooth closed hypersurface in 
IxE such that 

(i) the projections p : F — > X and 7r : F — > £ have rank n. This condition implies 
that the sets F (a) = ir^ 1 (a) , o G X and F (x) = p~ l (x) , x G X are hypersurfaces in 
X, respectively, in S. We call F incidence manifold. It can be defined locally by the 
equation / (x, a) = 0, where J is a smooth function in X x £ such that dJ ^ 0. If 
the hypersurface F is cooriented in X x S, a function / can be chosen globally; we call 
it incidence function. If F is not cooriented, one can choose local incidence functions 
I a ,a G A such that Iq = ±J a in the domain, where both functions /„, Ip are defined. 
The additional condition is 

(ii) the mapping q : F — » (X) is a local diffeomorphism, where G n_1 (X) = 
UxG" -1 and G™ -1 means the manifold of n — 1-subspaces in the tangent space T x of X 
at x; q (x, a) = (x, H) , where H denotes the tangent hyperplane to F (a) at x. It follows 
that for any point x G X and any tangent hyperplane H C T x there is locally only one 
hypersurface F (a) that contains x and is tangent to H. 



1 



Proposition 2.1 The conditions (i-ii) are equivalent to the inequality det <3> 7^ in F, 

where 
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where xi, ...,x n and <7i, are local coordinates in X and in £, respectively. 

Proof. Suppose that det $ 7^ 0. Then also d x I 7^ and d a I 7^ which implies (i). 
Choose a point (x , o"o) an d take a tangent vector 9 to F ((x ) at rr . The vector (6 1 , 0) is 
tangent to F at (xo,o- ) and the map q is well defined. Change the coordinates Xi and 
<jj,i = l,...,n in such a way that di (xo, o"o) jdo~i = di (xo, o"o) /dxi — 0,i — 2,...,n at 
this point. We have then 



det $ = 7— -— det^, * 



d 2 I 

dx i daj )ij , 



The inequality det^ 7^ implies that the forms dd x I (x , o) /da 2 , ...,dd x I (x ,a) /da n 
are independent. This means that the fields d/da 2 , ...,d/da n do not move the point x 
but rotate the tangent hyperplane to F (a ) at xq whereas the field d/da\ move the point 
xq. This yields (ii). The inverse statement can be proved on the same lines. ► 

It follows that the properties (i-ii) are symmetric with respect to X and S. The next 
condition is not symmetric: 

(iii) the projection p : F — ■> X is proper. If F satisfies (ii) and (iii) and the hyper- 
surface F (x) is not empty for a point x e X, then the mapping q : F (x) 
surjective, since the manifold G™ -1 is connected. 



G? 



is 



3 The Funk transform 

Consider a hypersurface F as above defined by a incidence function / that fulfils (i). 
Define the Funk (or Minkowski-Funk-Radon) transform for densities f in X with compact 
support by means of the integral 

s^O IE J\I(. :(T )\<£ J F (a) a xl 

where u — f / d x I is a n — 1-form such that d x I Aw = f. This form is defined up to a term 
d x I A Xi where x is a n ~ 2-form. Therefore the restriction of uo to the curve F {a) is a 
well-defined density. This density does not change, if we replace / by —1. Therefore the 
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Funk transform is well defined. The function Mi is also continuous. Suppose that the 
condition (iii) is fulfilled. If a density f is supported in a compact set K (s X, then Mi 
is supported in the compact set A = n (p^ 1 (K)) <e E. 

Example 1. Let X and E be unit spheres in Euclidean 3-spaces and the hypersurface 
F be defined by the global incidence function I (x, a) = X\<j\ + x 2 cr 2 + x 3 a 3 . The operator 
M coincides with the classical Minkowski- Funk transform, [3]. The Funk transform 
can be also defined on projective planes X = P 2 = S' 2 /Z2,E = P 2 , if we take ±7 as 
local incidence functions. The dual operator M° coincides with M through the natural 
isomorphism X = E. 

Example 2. Let (X,g) be a Riemannian 2- manifold with boundary and E be the 
family of closed geodesic curves 7. Take a density f = fdS, where dS is the Riemannian 
area form and / is a continuous function with compact support. Then we can write the 
geodesic integral transform as the Funk transform 

Mf( 7 )= f fds, 

if we take an incidence function I for the family F such that |V X /| = 1. Any smooth 
weight function w = w (x, a) can be included, by replacing the function I to w~ x F 

4 Above estimates 

The scale of Sobolev L2-norms ||-|| a ,a G M is defined for functions supported in an 
arbitrary compact set K C X. Fix a volume form dX in X and define ||f || a = ||/o|| 9a for a 
density f = f$dX with support in K. Denote by (X, fl) , (X) the space of densities 
(distributions), respectively, of (generalized) functions supported in K with finite norm 
|| - 1 1 . For an arbitrary compact set A C E we define the spaces H% (E, Q) , (E) in the 
same way. 

Proposition 4.1 For any family F that fulfils (i-ii), an arbitrary compact set K (s X , 
any smooth function e in E with compact support and any real a the inequality holds 

|| eMf ||* +(n-l)/2 < Ca || f ||a 

forieH%{x,n). 

Proof. The Funk transform can be expressed as an oscillatory integral 

Mi{a)= I exp(27m-I(x,o-))f (x) dr. 
Jk Jr 

The critical set of the phase function rl (x, a) is the hypersurface F (a) and the condition 
d x I 7^ implies that the phase function is non-degenerate. The corresponding conic 
Lagrange manifold is 

L = {(x, a, £, p) e T* (X x E) , / (x, a)=0,p = rd a I, £ = rdj, r ^ 0}. 
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Lemma 4.2 Rank of the matrix fj^y is equal to 2n in any point of A. 

Proof of Lemma. Suppose that the rank is less 2n. Then there exists a vector t = 
(Sx, 5a, 5£, 5p) in T* (X x E) tangent to L such that 5a = 0, 5p = 0. This yields 

dl (Sx) = 0, 5p = 5rd a I + rd x d a I (Sx) , 
8£ = SrdJ + rd 2 J (Sx) 

for a tangent vector St to BL The first line implies that the vector (tSx, St) fulfils the 
equation (tSx, St) $ = 0. By Proposition 12. II this vector vanishes, that is Sx = 0, St = 0. 
The second line gives <5£ = 0. ► 

By this Lemma the projections of L to T* (X) and to T* (E) are submersions. In other 
terms, L is locally the graph of a canonical transformation. The symbol a (x, a, £, p) — 1 
is a homogeneous function of £, p of order 0. The order m of the Fourier integral operator 
M satisfies the equation We have m + dimX x E/4 — N/2 = 0, where dimX x E = 2n 
and iV = 1 is the number of variables r. This yields m — (1 — n) /2, which means that 
the functional 



I (VO = / / / ex P (27m\Z~ (x, cr)) ^ (x, cr) dr 
is 7r 

defined for a smooth densities ip in X x E with compact support, is a distribution of the 
class j^ 1 -™)/ 2 (X x E,L) in the sense of Definition 25.4.9 of |S]. By Corollary 25.3.2 the 
operator eM defines a continuous map Hfr (X, O) -> Hl +{n ~ 1)/2 (E) for any real a, where 
A = supp e. ► 

Corollary 4.3 // F fulfils (i-iii), the Funk transform M can be extended to a bounded 
operator H« (X, fi) -> H* +(n ~ 1)/2 (E) /or any a, X and A = tt (p^ 1 (X)) . 

Due to (iii), we have supp Mf C A and the cutoff factor e can be dropped out. 



5 Dual Funk transform 

Let if be a density in E with compact support. Define the dual Funk transform as follows 

MV(x) = lim— / <p= I (2) 



e^O 2e J\I(x,-)\<e JF(x) &<tI 

If F satisfies (iii), then the manifold F (x) is compact for any x G X and the integral 
is well defined for any continuous density (p in E. The natural pairing 



(f,0)^(f,0)= f f0 
Jx 



is well defined for a densities f and functions on X provided one of them has compact 
support. 
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Proposition 5.1 The operator —M° is dual to M. 
Proof. We have 



JS </£ J F(a) ^-xl JF d x I Jf do-/ 

since dl = d x I + d a I = on F. The right-hand side equals 

" / f / -Tt=- I fM°(<p) = -(f,M°<p).+ 

JX JF{x) a <r J Jx 

6 Backprojection and two side estimates 

Definition. Fix some area forms dX in X and d£ in E. The back projection operator 

.\/ ' : // M (gdE) dX 

transforms functions defined in E to densities in X. 

Definition. We say that points x,y E X are conjugate with respect to F, if x ^ y 
and the form d a I (x, a) A d a I (y, a) defined in E vanishes. 

Theorem 6.1 If a family F has no conjugate points, fulfils (i-ii) and the condition: 
(iv) the projection q : tt" 1 (A) — > G n ~ l (K) is surjective for some sets K <<= X, A d E ; 
then for arbitrary cutoff function e such that e = 1 in A and any a > j3 the estimate 

||f f < C a \\eMif Hn - 1)/2 + Cp ||f f (3) 

holds for the Funk transform of densities f supported in K, where C a and Cp do not 
depend on f. 

Lemma 6.2 The composition M*eM is an elliptic PDO in K of order 1 — n. 

Proof of Lemma. Write f = fodX, where f is a function supported in K and calculate 

M - .Uf r e (a) dE r f(x) 

dX J F(y) d a I J F{a) d y I (x, a) 

= [ £{a)dJ: I dX /„(x) 

Jl(y,a)=0 d* 1 (V, <*) Jl(x,a)=0 d x J ( x i a ) 

e (a) dE 



h( y ,*)=o Ji(x,a)=o (r {y, o) A d a I {x, a) 
since dl = d x I + d a I = in F. We can write the right-hand side as J A (y, x) f (x) , where 

A(yx) = -[ " dS (4) 

JF(x)nF(y) d CT I (y, a) A dj (x, a) 
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The dominator does not vanishes for y ^ x, since there is no conjugate points. The 
quotient Q in (@J is well defined as n — 2-form up to an additive term d a I (y, a) A 
S + d a I (x, a) A R, where S and R are some n — 3-forms. The integral of this term 
along the smooth manifold F (x) fl F (y) vanishes and the function A (y, x) is a well 
defined and smooth, except for the diagonal. Near the diagonal we can write / (y, a) = 
I (x,a) +J2(Vi- x i) 91 (x,a) /dxi + O (\y - x\ 2 ) and 

dj (x, a) A dj (y, a) = ^ (Vi ~ x i) ( x > a ) A d ° — ^7 VO[\y- x\ 2 ) . 

i % 

The forms d a I (x, a)Ad a dI (x, a) /dxi, i = 1, n do not vanish and are linearly indepen- 
dent, since of Proposition 12.11 Therefore the product d a I (x, a) A d a I (y, a) is bounded 
by c \x — y \ from below as y — > x. Therefore we have A (y, x) = a (y) \x — y\~ +0 (1) near 
the diagonal, where a is a smooth positive function. This implies that M*sM is a classical 
integral operator on K with weak singularity, moreover it is a pseudodifferential operator 
of order 1 — n. It is an elliptic operator, since the symbol a is positive. ► 

Proof of Theorem. The support of the function M*eMi is contained in the compact 
set p (suppe) C X. By Proposition 14. II M* is (n — 1) /2-smoothing operator, which yields 

\\M*eMf ir +n_1 < C ||eMf |r +(n ^ 1)/2 . (5) 

By Lemma 16.21 the operator M*eM is elliptic of order n — 1, therefore the standard 
inequality holds 

||f || a < C a \\M*sMi\\ a+n - 1 + C p ||f f 
for an arbitrary f3 and some constants C a ,Cp. Taking in account (0) yields (jSJ). ► 

Corollary 6.3 The eigenvalues of the operator M*eM numbered in decreasing order 
satisfy the estimate 

cjfc(l-»)/2 < \ k < Ck^l\ k > k . 

For the Radon transform the eigenvalues are calculated in [TU] . 

Corollary 6.4 Suppose that for some (3 < a the equation Mi = 0, f € (X, Q) implies 
f = 0. Then the two-side estimate holds: 

c a iif ir < \\eMif +{n - i)/2 < c a iifir . (6) 

Proof. The right-hand side inequality follows from Proposition 14.11 Suppose that the 
left-hand side estimate does hold for no c a . Then there exists a sequence {4} C Fl% (X, Q) 
such that 

HffclL > k \\Mi N \\ a+l/2 , ||f fc || ^ = 1, k = 1, 2, ... (7) 
The inequality ([3]) implies that ||ffc|| Q < 2Cp for k > 2C a and ||Mffc|| a+1 / 2 — > 0. Because 
the imbedding (X, Vt) — > H^ K (X, Q) is compact, we can choose a subsequence (denote 
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it again {ffc}) such that — > g in (X, ft) . By Proposition ^. II || ML — > ^g|| /3+ ( n _i)/ 2 — » 
0, which implies Mg = 0. By the condition g = 0; it follows that ||ffe|L —> in contradic- 
tion with Q. ► 

Remarks. Mukhometov's result [9] implies the estimate ||f ||° < C \\Mi\\ l for the case 
n — 2. An estimate of this kind for more general situation was obtained by Sharafutdiniv 
|15j . Ch. IV. Inequalities for Sobolev norms are well known for the Radon transform. 
Estimates for shift derivatives of order a + 1/2 (n = 2) were obtained by several authors. 
Natterer [10] has shown that (EI) holds also for angular derivatives. For the attenuated 
Radon transform see Rullgard |14j . 

Our approach is similar to that of Lavrent'ev and Bukhgeim [7], where the composition 
M*M was described as an integral operator in the local case. Guillemin [4], [5] has defined 
the 'generalized Radon transform' R for an arbitrary double fibration. This transfrom is 
treated as an elliptic Fourier integral operator and R*R is shown to be a pseudo differential 
elliptic operator under the 'Bolker condition'. This condition is equivalent to absence of 
conjugate points in our situation. More details are given in the paper of T. Quinto [T3] . 



7 Range conditions and approximation 

Let K be a compact set in X and a G M. We define H a (K, ft) to be the dual space 
of H^ a (X) and use the notation ||-|| a for the norm in H a (X, ft) . The trace operator 
HI (X, ft) — > H a (K, ft) is well defined and is open for an arbitrary compact set L C X 
such that K d L, since H^ a (X) is a subspace of H£ a {X) . Therefore H a (K, Q) can be 
realized as the quotient space of H£ (X, Q) modulo the kernel of the trace operator. The 
last one consists of densities / supported in L\K. For any (3 > a we have the operator 
r]' : H 13 (K, Cl) — > H a (K, Q,) , which is dual to the natural imbedding i] : H^ a (X) — > 
H^f (X) . If the boundary of K is smooth, the imbedding r] has dense image and 7/ is 
injective. Then we can define the intersection if 00 (K, Q) = T\ a H a (K, Q) ; any density f G 
H°° (K, O) is smooth in the interior of K. Similarly, we define H a (K) = (H^ (X, fl))' . 

Suppose that the incidence manifold F fulfils (i-iii). If a density f is supported in a 
compact set K, then the support of Mf is contained in the compact set A = tc (p -1 (K)) . 
The hypersurface F (x) is compact for any point x G X and the dual transform M° is well 
defined for all continuous densities in S. Moreover, it can be extended to a continuous 
operator M° : H a (A, Q) -> ff a +(»- 1 )/ 2 (K) for any a by means of the duality 

(M°g, f ) = - (g, Mi) , f G H K a - {n - 1)/2 (X, n),gEH a (A, Q) . 

By Proposition 14.11 we have Mf G H^ a (E) , hence the right-hand side is well defined. If 
a > 0, the function M°g defined by this formula is equal to the integral (|2J), which has 
sense, at least, for almost all x G K. 

Theorem 7.1 Suppose that F satisfies (i-iii) and has no conjugate points. Then for 
any K d= X and arbitrary a G IRU {00} t/ie image of the Funk operator 

M : H« (X, ft) -> H a h Hn ' 1)l2 (S) , A = 7T (p" 1 (X)) 



7 



is closed and coincides with the subspace of functions <p £ if^ +(n (£) swc/i i/iai J gip = 
/or any solution g £ H~ a ~^ n ~ 1 ^ 2 (A, fl) 0/ t/ie equation 

M°g(x) = 0,x £ K (8) 



Proof. The image of M is closed by Theorem 16.11 thereby it coincides with the polar 
of the kernel of the dual operator M°. ► 

Theorem 7.2 If F fulfils (i-iii) and has no conjugate points. Then for any set K <s 
X with smooth boundary and arbitrary real a £ M any density g £ H a (A, Q) , A = 
7r (p^ 1 (K)) that fulfils (0§P can be approximated by solutions h £ H°° (A, Q) = r\ a H a (A, fi) . 

Proof. Let Sol 13 denote the space of solutions of (|S]) in the class (A, fi) . We show 
first that g can be approximated by elements of Sol' 3 for any f3 > a. It is sufficient to 
check that any functional on H a (A, fi) that is equal to zero on Sol 5 also vanishes on 
g. The dual space is isomorphic to H^ a (£) , which implies <ft £ ^a" (^) • By Corollary 
14.31 the Funk transform defines the continuous operator 

M p :H-/- {n - l)/2 (X,fi)-ff/(E). 

By Theorem 17.14 the range of this operator is closed and coincides with the polar set of 
Sof 3 . It follows that, (j) = Mip for a density £ H~ P ~ {n ~ 1)/2 ^ ^ _ By Theorem Owe 
have V e tf^"- 1 )/ 2 (X, ft) and can write 

(0,g) = (M^,g) = -(^M o g) = O, 

since the function M°g is well defined as element of jj^^ 1 ^ 2 (X) . This yields that g 
is contained in the closure of the space Sol . 

Now we approximate g by elements of the space H°° (A, Q) . Let ||-||^ be the norm 
in the space H a (A, Q) , which is dual to the norm We may assume that the 

norm ||-||^ is monotone increasing function of a. Take an arbitrary e and choose a 
function hi £ H a+1 (A, Q) such ||hi — g||^ < e/2, then we choose a function h 2 £ 
ffa+2 (A, fi) such that ||h 2 — hi||^ +1 < e/4 and so on. We obtain a sequence {h^} such 
that ||hfc + i — hfc||^ +A; < 2~ k ~ 1 e for k = 1,2, ... This sequence converges to an element h in 
any space H p (A, fi) , f3 > a. It follows that h £ H°° (A, Q) and ||h - g||° < e. ► 



Corollary 7.3 Under conditions of Theorem \ 7. 1\ it is sufficient to check the equation 



J gip = for densities g £ H°° (A, Q) that satisfies 



8 Thermoacoustic tomography 

We apply the above results for the thermo/opto/photoacoustic geometry. First, consider 
the case of complete acquisition geometry. Let X be the open unit ball in an Euclidean 
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space E, Sr be the sphere of radius R > 1 and E = Sr x JR.. The manifold FcIxE 
is given by the equation I(x;y,r) = \y — x\ —r = 0,y £ Sr,0 < r. The manifold F 
obviously fulfils (i),(iii) and has no conjugate points. Check that the condition (ii) is 
also satisfied. It is sufficient to check that, the sphere F (y, r) can not be tangent to 
F(z,s) at a point x £ X, if the points (y,r) and (z, s) are sufficiently close in E. The 
condition \y + z| > 2 is sufficient for this. The Funk operator 

Mi(y,r)= I f (x)dS, f = f dx (9) 

J\x-y\=r 

is the spherical integral transform, where dS* is the Euclidean surface area form on spheres. 
The kernel of dual transform M° consists of densities ip = (pdSdr in E such that 

J F(x) d a I J d(\x-y\-r) J F{x) 
for any x £ X. Theorem 17.21 yields 

Corollary 8.1 For any compact set K C X with smooth boundary and arbitrary a £ R 
U {oo} the range of the Funk operator M : (X, Q) — > H^ + ^ n (E) coincides with 
the set of functions g in E such that 

fw = o (10) 

for any density ip = ^dS'dr such that <p £ C°° (E) and 

[ <f)dS = 0, x £ K. (11) 

Jf{x) 

Remark. For the operator M acting on C°°-densities the range conditions were given 
in the papers [11] , [1] , [2] . The conditions of pQ and |2] give full description of the range 
of M, but have implicit form. 

We extract some explicit range conditions from Corollary 18. II For an arbitrary x £ K, 
the manifold F (x) is the intersection of the cone surface \y — x\ = r with the cylinder E. 
This intersection is contained in the hyperplane 

P ( x ) = {y, s; 2(x,y) + s = \x\ 2 + R 2 } cExK, 

where we set s = r 2 . Thus, the condition (fTTI) means vanishing of integrals of <pdS over 
intersections of E with the hyperplanes P (x) , x £ K. Suppose that is a polynomial in 

s : (y, s) = 4>k (y) (s - R 2 ) k and have 

/ <J>dS = J2 I <f> k (y)(\x\ 2 -2(x,y)) k dS(y). 

JF(x) £ Js 
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Set x = tz for \z\ — 1 and < t < 1 and develop the right-hand side in powers of t 



[ 0dS = V / <j> k (y)(t 2 -2t(z,y)) k dS 

JF{x) Js 

r <f> (y)dS-t [ 2<f> 1 (y)(z,y)dS 



+ t 2 j [h(y) + ^2(y)(z,y} 2 ]dS 

-t 3 J [40 2 (y) (z, y) + 80s (y) (z, y) 3 } dS 

+ 1 4 J [02 (y) + 1203 (y) (z, y) 2 + 160 4 (y) (z, yf] dS 
+.. =0 

The right-hand side vanishes for all t, which yields the system of equations 

0o (y)dS' = 0, 
0! (y) (z,y)dS = 0, 
[402 (y) (z,y) 2 + <p 1 (y)] dS = 0, 



[203 (y) {z,yY + <t> 2 (y) (z,y)] dS = 0, 
[160 4 (y) (z, yf + 1203 (y) (z, y) 2 + 2 (y)} dS = 0, 



Corollary 8.2 Any solution (0o, 0i, 02, •••) of this system such that <pj = for all j > k 
for some k yields a function (y, s) = ^ <pj (y) (s — R 2 ) 3 that is a polynomial in s of 
order k, fulfils ( fill) and is orthogonal to the range of M. 



There are many solutions of this form, since the system has triangle form with diagonal 
terms 

0j (y) (z,y) J dS, \z\ = 1, j = 0,1,..., k. 



To solve these equation we only need to fix the moments of 0j of degree j. There are only 
( n+ ^7 ) nnearr y independent j-moments, hence one can find infinitely many independent 
solutions which are finite sums of harmonics. In particular, we can take for O any function 
on the sphere with zero average, an arbitrary function <pi with zero linear moments of 
0i and set <pk = for k > 1 etc. The range conditions of S. Patch [11] are apparently 
contained in (TTOT) for polynomial 0. 
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9 Partial scan and Kaczmarz method 



In the case of the partial scan geometry the analysis is more complicated. Let again X 
be the open ball of radius 1 in E and £5 = S$ x R +! where S$ = {y; \y\ — R, y\ > —6} 
for some 5 > 0. The manifold F$ is defined in X x £5 by the same incidence function 
J, that is, Fg(y,r) is the sphere of radius r with the center y G S5. The manifold F5 
fulfils (i),(ii) and has no conjugate points, but does not fulfil the condition (iii). On 
the other hand, F$ satisfies (iv) for the unit half-ball K = {x, \x\ < l,x± > 0} and A = 
{y, r : R— I < r < R + 1}. Consider the Funk transform 

M : H« (X, ft) -> Hl +{n ' l)/a (£) (12) 

defined as in ([9]). Take a smooth function Eq > on R supported in the interval (—5, 00) 
such that £0 (£) = 1 for < t < 1. Set e (y) = e (2/1) and consider the operator M*eM. 

Proposition 9.1 T/ie operator [TB) is injective for any 5 > and arbitrary a > 1/2. 
TTie inequality holds for the operator M*eM and any a > 1/2. 

A We prove that the equation Mi = in S e for a density f G (X, Q) implies 
f 7^ 0. This condition means that the spherical means of f vanish for spheres centered 
at points y G S$. By the Lin-Pinkus theorem [8] this implies that either f = or a non- 
trivial harmonic polynomial h vanishes on Ss (the continuity condition for f in [8] can 
be weakened). The last option is impossible, since Ss is strictly convex and Ah can not 
vanish identically near Sg- The estimate ([6]) now follows from Theorem 16. II and Corollary 

ES ► 

We show that the Kaczmarz method can be adapted for inversion of the operator ffl2|) 
as well as for a general Funk operator M. Let K, A be compact manifolds, occasionally 
with boundaries, F be a closed hypersurface inKxA that fulfils (i) and (ii). We want 
to find a solution f G H° (X, fl) of the equation 

Mi = cp (13) 

for a function ip G H°(A). By Proposition 14.11 the left-hand side is contained in the space 

R (n-l)/2 ( A ) c H (A) . 

Example 3. Take for K the closed unit ball in E, for A the product Ss x [R — 1, R + 1] 
and for i 7 " the manifold of spheres as above. The conditions (iO and (ii) are fulfilled. 

Fix a volume form dX in K, a volume form dS in A and consider the operator 
MM* : H° (A) -> #° (A) . It is non-positive; set R = -MM* + 61, where / is the 
identity operator and 9 > 0. The operator R is self-adjoint, positive and invertible. 
Following [TU], we choose a real parameter u; and set Q = I — uoM*R~ 1 M. We use the 
notation ||-|| = ||-||° . 

Lemma 9.2 We have \\Qg\\ < ||g|| for < to < 2 and any g G H° (K,Q) such that 
Mg^O 



11 



Proof. We have by Proposition 15.11 

||Qg|| 2 = || g || 2 _ 2u (g, M*R~ 1 Mg) + oo 2 (M*R~ 1 Mg, M*R~ l Mg) 
= ||g|| 2 - 2u (Mg, R~ x Mg) - u 2 (R^Mg, M*MR~ 1 Mg) 
= ||g|| 2 - 2oo (Mg, R^Mg) + u 2 (R' l Mg, Mg) - u 2 (R~ l Mg, R~ x Mg) 
< \\g\\ 2 -uo{2-uo)(Mg,R- 1 Mg), 

since (R^Mg, R~ l Mg) > 0. The term (Mg, R~ l Mg) is positive, if Mg ^ 0. ► 

Take an arbitrary density f° and construct the sequence f fc , k = 1, 2, ... by means of 
the recurrent formula 

i k + l = f fc + uM*R- 1 (ip - Mi k ) . 

Theorem 9.3 If M is injective and (p fulfils the range conditions, we have i k —>■ i, where 
f is a solution of (TT3|). 

Proof. We have 

Q (f k - f) = f k - ooM*R~ 1 Mf k - f + uM*R~ 1 M{ 

= i k + uM*R- 1 ((p - Mf k ) - f = f k+1 - f . 

It follows that 

||f fc+1 -f|| < ||f fc -f|| < ... < ||f°-f|| 

and f k — > g strongly in L 2 (K, Q) . We have \\Q (g — f) || = ||g — f || , which yields g = f by 
Lemma [9.21 ► 

Another inversion method is developed by Popov and Sushko [12] . 
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